Ultracold Bosons in a Tilted Multi-level Double- Well Potential 
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The iV-body problem in a double well requires new features for quantum information processing, 
macroscopic quantum superposition, and other fundamental studies of quantum many body physics 
in ultracold atoms. One needs (a) tilt, and (b) to go beyond the single-particle ground state in each 
well, i.e., to two or more energy levels. For (a), we show that a small potential difference between the 
wells, or tilt, causes the decoherence of Schrodinger cat states. However, these states reappear when 
the tilt can be compensated by atom-atom interactions; these tilted cat states constitute partial cats 
that are protected from decoherence by the many body wavefunction. For (b), we provide explicit 
criteria for when two energy levels are needed to describe the state space. For typical experimental 
parameters, two levels are indeed required for creation of cat states. 

PACS numbers: 



Recently, Bosc-Einstcin condensates (BECs) in double- 
well potentials have been the subject of diverse and ex- 
citing research. For instance, such systems can be used 
to search for deviations from Newtonian gravity at small 
distances [l|, Q and to store and retrieve optical informa- 
tion 0] . Interest in these systems is not limited to practi- 
cal applications; BECs in double-well potentials also pro- 
vide an ideal medium for the study of fundamental quan- 
tum many-body phenomena, such as macroscopic quan- 
tum tunneling 0,11, H0B (3 an d macroscopic superpo- 



sition states, or "Schrodinger cat" states 0, [ljj 11, 
To describe these double-well systems, a Hubbard-like 
Hamiltonian has often been employed [HI, 0, EH ■ How- 
ever, these studies not only assume that the trapping po- 
tential is symmetric, but effects of excited levels are com- 
pletely neglected. Variational methods have indicated 
that excited levels play a significant role 
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Moreover, a multi-level picture of a tilted double- well po- 
tential is required for the creation of a quantum computer 
from neutral atoms (l9l . |20| , atom-chip-based gravity sen- 
sors [l], [H, and in_the study of quantum transport phe- 




231 ] - Thus, the current experimental con- 
text of double- well potentials has created an urgent need 
for a new theoretical analysis of the many-body double- 
well problem. 
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FIG. 1: (color online) Energy eigenstates of a symmetric 
double well. Shown are the probability amplitudes |c£ | for 
a small range of Fock index n and energy eigenstate index k. 
Amplitude is indicated by hue. The high barrier limit in which 
macroscopic superposition states dominate the spectrum is 
depicted. The total number of atoms is N = 20; above the 
dotted line at Fock index n = 20, one atom occupies the 
higher energy level; above the dotted line at n = 60 two 
atoms occupy the higher level; etc. Inset: schematic of the 
double-well potential. 



In this Letter we use a two-level Hamiltonian to in- 
vestigate the stationary states of a BEC in a tilted, one- 
dimensional (ID), double- well potential. A schematic of 
our potential is shown in the inset of Fig. [1] In addition 
to interaction with the electromagnetic vacuum 01 , ther- 
mal effects 
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I , and dissipation [121 ] , we find that tilt also 
causes the collapse of cat states when the barrier is high, 
an effect we term potential decoherence. Unlike for other 
forms of decoherence, such states reappear when the tilt 
is compensated by atom-atom interactions. We call this 
a tunneling resonance because tilt suppresses tunneling 
between wells except when these states reappear 24 , 25| . 
Finally, we present novel, formal bounds for the use of 
a one-level approximation. This approximation is typi- 



cally thought to hold when the interaction energy is much 
smaller than the energy level difference [l8||. However, we 
show that the effects of the excited energy level cannot 
be neglected even in this regime. 

Past studies of the iV-body double-well problem have 
focused on symmetric traps with one allowed energy level. 
Such systems map onto other physical problems, includ- 
ing a classical nonrigid pendulum Q and the Lipkin- 
Mcshkov-Glick model from nuclear physics [26|. More- 
over, a double- well potential coupled to a heat bath can 
be used to study vibration modes in water molecules [27] ] . 
For the symmetric one-level model, tunneling dynam- 
ics [H, 0, El, HH and the response to changes in experi- 
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mentalparameters such as barrier height have been ana- 
lyzed [Hi 15, [2^| in various regimes. The behavior 



of the ground state energy in a tilted double-well has 
been studied in the one- level approximation [24|; how- 
ever, applications such as quantum information process- 
2C 
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call for a two-level description. In our study of 
the double- well potential, we relax two assumptions com- 
monly made in previous studies of similar systems: the 
symmetric trap assumption and the one-level assump- 
tion. This leads to two new energy scales in our problem. 

An approximate two-level Hamiltonian for N weakly 
interacting bosons in a tilted double-well potential is 



H = H° + H 1 +H ( 



(1) 



where 



-(AV/2) (ni-n^+E^ni 



is the usual one-level Hamiltonian and 



H m = U m Y. 



(2) 



(3) 



couples the energy levels. Here [S^,^-/ ] 
= [UX]=0,aDdrt 
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- bjbj. Equation fTJ) can 
be derived from first principles quantum field theory for 
weakly interacting bosons at zero temperature [29j . The 
superscripts £, £' € {0, 1} are the energy level indices, the 
subscripts j,j' € {L,R} are the well or site indices, J e 
are the tunneling energies, U e and U 01 are the interac- 
tion energies, E e is the energy of the ^th excited level, 
and AV is the tilt. Setting E° = 0, the energy difference 
between levels, or level spacing, is E 1 = Two. The exten- 
sion of Eq. ([T]) to an infinite number of sites leads to the 
two-band Bosc-Hubbard Hamiltonian. 

The two-level Hamiltonian allows for on-site interac- 
tions, tunneling between wells, and hopping between lev- 
els. The tunneling terms J 1 allow single particles to 
tunnel between wells in the same level. Pairs of par- 
ticles in the same well interact with interaction energy 
U e if they are in the same energy level and U 01 if they 
arc in different levels. Interactions can be either repul- 
sive, U l > 0, or attractive, U £ < 0. Furthermore, while 
single-atom transitions between energy levels are forbid- 
den, two atoms can hop together between energy levels 
with amplitude U 01 . Thus the energy levels are cou- 
pled by the inter- level interaction energy U . Interac- 
tions between atoms in different wells are much smaller 
than U° and have been neglected in Eq. (fT]). The pa- 
rameters J £ , fku, and U are determined by overlap inte- 
grals of the localized single-particle wavcf unctions. Left- 
and right-localized wavefunctions are constructed by su- 
perpositions of the appropriate symmetric and antisym- 
metric eigenfunctions of the single-particle Hamiltonian. 



While the interaction and tunneling energies are inde- 
pendent parameters, J e and fku arc not. The tunnel- 
ing and interaction energies satisfy J° -C J 1 <C hu, 
U 1 = (3/4)[/° and U 01 = (1/2)17°. In Figs. 1, 2, and 
3 we set J°/hu> = 4 x 1(T 7 and J 1 /Hlu = 3 x 1(T 5 . 
Finally, we note that the use of single-particle wavefunc- 
tions and only a few levels is appropriate to the regime 
\NU°\ < 2fko. When \NU°\ > 2fko our approximation 
is inaccurate and alternative treatments become neces- 
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An arbitrary state vector in Fock space is given by 



I*) 



En= c n \n), |n)e{K,0®K,n^)}, (4) 



where n £ is the number of particles in the ^th level of the 
jth well. We require the total number of particles in the 



rij, to be constant. The Fock in- 
„o , „i 



double-well, N 

dex n increases with n^+n\, and n L . The dimension 
of the Hilbert space is Q = (N + 3) (A + 2)(N + l)/6; in 
the one-level approximation this reduces to simply 2V+1. 
The energy eigenstates and eigenvalues of the Hamilto- 
nian (fT]) are given by H\4>k) = £fc|</>fc), for < A; < fi — 1. 
The probability amplitudes are defined as cl = (n\<j)k} 
for < n < — 1. Throughout our treatment, we con- 
sider the regime |JVC/°| <?C 2frui. Furthermore, we work in 
the high barrier limit J° -C | f7°| , as it is instrumental to 
the creation of stationary cat states. For simplicity, we 
restrict our discussion to repulsive interactions, U° > 
and positive tilt, AV > 0. However, our results hold for 
U° < and AV < as well. 

Characteristic eigenstate amplitudes arc shown in 
Fig.[JJfor a symmetric trap. For a given n, the number of 
particles in the excited level is equal to the number of hor- 
izontal dashed lines below n. In the case of Fig.[TJ mixing 
of levels is negligible. The lowest excited group of eigen- 
states, for instance, corresponds to full occupation of the 
lowest level. Likewise, the eigenvalues occur in N + 1 
groups approximately separated by the level spacing Hlj. 
As the number of particles, the interaction energy, or the 
tilt increases, the spacing between groups decreases, as 
in Figs. HJa) and^b). In a symmetric trap, the first of 
many eigenvalue crossings occurs when the condition 



^ ^crit 



2Huj/(N 2 - 1), 



(5) 



is violated. When U° > U® Tit , cat states with nonzero 
occupation of the excited level emerge among the first 
A + 1 eigenstates and a two- level approximation must be 
used. As the interactions increase, the lowest eigenstates 
include cat states with successively larger occupation of 
the excited level. Eventually, even the ground state will 
have significant contributions from the excited level. In 
a tilted potential, condition ([5]) becomes 



AV<AV c , Ti t = 



'tku/N, 



2U° 



y/l + 2{huj/U )-N 



u°<u° 
u°>u° 



(6) 
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FIG. 2: Breakdown of the one-level approximation. Energy 
eigenvalues versus (a) interaction energy with AV = and 
(b) tilt with J°/U° = 0.1 for a ten-atom system. The dashed 
red lines in panels (a) and (b) correspond to U° = U® lit and 
AV = AVorit, respectively. Eigenvalue crossings cause states 
with nonzero occupation of the excited energy level to emerge 
among the N + 1 lowest-lying eigenstates. 



FIG. 3: Potential decoherence and tunneling resonances. 
Eigenstate probability amplitudes for a small range of Fock in- 
dex n and eigenstate index k for N = 10, J°/U° = 0.1, and (a) 
AV/U° = 0, (b) AV/U° = 10~ 2 , and (c) AV/U° = 6. (b) A 
small tilt collapses the stationary cat states, (d) Avoided 
crossings in the energy eigenvalues indicate a reappearance of 
partial cat states, as in (c). 



where we have defined t/° ax = U° rit (N +1)/(AN). These 
results can be obtained formally in the limit of small tun- 
neling. In the opposite, noninteracting limit, we have 
developed a formal criterion involving the hopping as 



well [31|: N < 1/2+ (Hlj - J 1 )/(2J°). 

While direct product states are not true eigenstates of 
the two- level Hamiltonian for U° > 0, the N + 1 lowest 
excited eigenstates have negligible contributions from the 
excited level when conditions ([5]) and © are met. In a 
symmetric potential, these states are of the form 

^±;«/> = l*±;v>®|o,o) ) (7) 

to (N-2v- l)th order in J°/U° for < v < N/2. Here 



|*±; v) = (h N - v) ± \N - v, v))/y/2, 



(8) 



represent partial cat states. To (N — 2u)th order in 
J° /U°, the energy difference Ae„ of antisymmetric (-) 
and symmetric (+) pairs of states is 



Ae u = 



4U Q [J /{2U°)} N - 2 »(N -v) 
v\[{N -2v- l)!] 2 



(9) 



which is a very small number. In the presence of a high 
barrier, the eigenstates therefore occur in nearly degen- 
erate pairs of entangled states 0, [T^j when the poten- 
tial is symmetric. Here v = represents the extreme 
cat state in which all atoms simultaneously occupy both 
wells. Characteristic probability amplitudes for the 7V + 1 
lowest eigenstates are shown in Fig. EJa). 



Because the level splitting between antisymmetric and 
symmetric pairs is so small, small perturbations can mix 
these states 0] and produce a localized state of the form 
1 1/, N — v) (g) |0, 0). Indeed, cat states are highly sensitive 
to tilt AV, as illustrated in Fig. [31(b) . The cat states of 
Eq. (J7J) are destroyed when 



AV > 2Ae u / (N - 2v) 



(10) 



Small imperfections in the external potential thus con- 
stitute a source of quantum decoherence. In addition 
to dissipation and measurement, potential decoherence 
therefore poses a further difficulty in the engineering of 
cat states in experiments. Because condition (|10p is min- 
imized when v = 0, extreme cat states are the most sen- 
sitive to imperfections in the double-well, making them 
an unlikely candidate for experiments. Partial cat states, 
v > 0, on the other hand, are more robust with respect 
to potential decoherence [32j. 

Despite their fragility, cat states reappear periodically 
for certain values of the tilt. Such resonances occur when 

AV = AV P = 2 P U°, p e {1, 2, . . . , N - 1}. (11) 

In this case, the potential difference can be exactly com- 
pensated by the interaction of p atoms in the lower well. 
To (TV — 2v — p — l)th order in J° /U , the energy eigen- 
states are partial cat states of the form 

\<f>±;v,p) = \^±\v,p}(»\0, 0), (12) 
\^ ± ;p,p) = (\p, N-v)±\N-v-p, v+p))/V2, (13) 
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for < v < (N — p)/2. The reappearance of cat states 
is shown in Fig.[3jc) for AV = AV3, i.e., the third reso- 
nance. Because these states also occur in nearly degen- 
erate pairs, the tunneling resonances are easily identified 
by avoided crossings in the energy eigenvalues, such as 
those displayed in Fig. (2Ud). At odd integer values of 
AV/U°, the eigenstates are maximally localized with the 
largest energy splitting between pairs of states. Near a 
resonance, the eigenstates become localized when 



\AV- AV P \ >2AeP/{N-2is- P ), 



(14) 



where As? oc (J° /U°) N ~ 2 "~ P is the energy difference be- 
tween the states |0_; v, p) and \<j> + ; v, p). For the special 
case u = 0, we find 



Ael 



AU°[J°/{2U°)] n -p(N -p) 
(JV-p- 1)! 



JV! 



p\(N-p)V 



(15) 



to (JV — p)th order in J°/U°. The understanding of tun- 
neling resonances is vital in systems in which tilt is ap- 
plied deliberately. 

To demonstrate the robustness of the resonant cat 
states, we consider JV = 100 87 Rb atoms in a ID ana- 
log of the double- well potential of Ref. [3] : 

V(x) = -V! cos 2 (2kx) - v 2 cos 4 (fca; - tt/4 - 9), (16) 

for kx £ [-7r/4,37r/4], v 2 / Vl < 2, and 9 € [0,tt/4]. 
Here k = 2tt/X where A = 810 nm. We set v\ = 
V2/0.15 = 106£ r where E r /h = 22 kHz is the recoil 
energy of the lattice (33|. The radial trapping frequency 
is u>± = 3.2 kHz. For such a potential, fuj/E r = 36, 
NU a /(2huj) = 0.10, and J°/U a = 0.10. The criti- 
cal interaction energy and tilt are U° Tit /U° = 0.10 and 
AV cy ;it/E r = 0.36, respectively. In a symmetric trap, 
9 = 0, the extreme cat states \(j>±;0) will collapse for de- 
viations in the tilt on the order of lO -287 ^,.. However, at 
the 98th resonance, i.e., when AV = AVgg = l^E r , the 
few-atom superposition state |<J>±;0, 98) can withstand 
deviations in the tilt up to ±0.098i5 r . The 98th reso- 
nance corresponds to 9 = 0.52 ± 1.2% (34|. Clearly, the 
many-body wavefunction protects partial cat states; a su- 
perposition state of a few atoms without the presence of 
a many-body cushion is not experimentally realistic [25| . 

In conclusion, we have shown that not one, but three 
energy scales are required to describe macroscopic su- 
perposition states of ultracold atoms in a double well: 
hopping, tilt, and energy level spacing, all in ratio to in- 
teraction. These energy scales arise naturally within the 
context of quantum information processing with neutral 
atoms, where they are required for quantum operations. 
After solving the appropriate Hamiltonian, we showed 
that tilt is a prevalent source of quantum decoherence 
which cannot be neglected in experiments. Moreover, we 
provided simple, formally derived criteria for when all 
three energy scales are required. Our treatment of the 



two-level double-well potential was restricted to ID for 
simplicity. While all considerations made here also apply 
to 2D and 3D, there is one very important difference. In 
the latter case, angular momentum associated with the 
excited state in each well introduces another quantum 
number, as we analyze in detail elsewhere 3l|. 
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